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This paper briefly summarizes the search for astronomical sources of high-energy neutrinos using the AMANDA- 
BIO detector. The complete data set from 1997 was analyzed. For _E,i > 10 TeV, the detector exceeds 10,000 
m^ in effective area between declinations of 25 and 90 degrees. Neutrinos generated in the atmosphere by cosmic 
ray interactions were used to verify the overall sensitivity of the detector. The absolute pointing accuracy and 
angular resolution has been confirmed by the analysis of coincident events between the SPASE air shower array 
and the AMANDA detector. Preliminary flux limits from point source candidates are presented. For declinations 
larger than -1-45 degrees, our results compare favorably to existing limits for sources in the Southern sky. We also 
present the current status of the searches for high energy neutrino emission from diffusely distributed sources, 
GRBs, and WIMPs from the center of the earth. 



1. INTRODUCTION 

The AMANDA-BIO high energy neutrino de- 
tector was constructed between 1500-2000 m be- 
low the surface of the Antarctic ice sheet where 
the optical properties are suitable for track re- 
construction R]. The instrumented volume forms 
a cylinder with outer diameter of 120 m. The sur- 
face electronics are located within a kilometer of 
the Amundsen-Scott Research Station at the ge- 
ographic south pole. The detector was commis- 



sioned in February lOOTMH], and initial scientific 
results were presented at the XXIVth Interna- 
tional Cosmic Ray Conference H. Reconstruction 
methods and detector calibration techniques are 
introduced in reference [g. 

AMANDA-BIO consists of 302 optical modules 
(OMs) that contain an 8 inch diameter photomul- 
tiplier tube controlled by passive electronics and 
housed in a glass pressure vessel. They are con- 
nected to the surface by an electrical cable that 



provides high voltage and transmits the signals 
from the OM. The simple, reliable system archi- 
tecture is responsible for the low fraction of OM 
failure (< 10% after several years of operation, 
although most of the failures occur within a week 
of deployment). 

In January, 2000, AMANDA-II was completed. 
It consists of 19 strings with at total of 677 OMs 
arranged in concentric circles, with the ten strings 
from AMANDA-BIO forming the central core of 
the new detector. New surface electronics consol- 
idates several triggering functions and adds func- 
tionality. New scalers were installed that provides 
millisecond resolution - important for Supernova 
studies. Several technologies were deployed to 
evaluate their utility and readiness for future ex- 
pansion to larger systems. 

The search for astrophysical sources of high en- 
ergy {El, > 1 TeV) neutrinos is one of the cen- 
tral missions of the AMANDA detector. Figure |l| 
provides a survey of model predictions for the 
flux of high energy neutrinos from point objects. 
It also contains the flux limit reported in this 
manuscript for sky bins with declinations greater 
than 30 degrees. Most theoretical models of po- 
tential astrophysical sources of neutrinos predict 
that the energy spectrum is very hard, approxi- 
mately E~^ . Due to the hard energy spectrum, 
the most probable energy of the detected neutrino 
is well above 1 TeV (typically 10-30 TeV). 

In this report, we describe the status of a gen- 
eral search for continuous emission from point 
sources in the northern sky, restricted to declina- 
tion greater than -1-5°. The observation of atmo- 
spheric neutrinos and detailed studies of downgo- 
ing atmospheric muons have been used to study 
the sensitivity of the detector. The observed dif- 
ferences between observation and prediction are 
within ±30%, which is within the systematic un- 
certainty. By restricting the atmospheric neu- 
trino analysis to a search cone with a half angle of 
15° about zenith, a search was made for neutrino 
emission from WIMP annihilation near the center 
of earth. Finally, we report preliminary results for 
a search for high energy neutrinos from Gamma 
Ray Bursts (GRBs). Neutrino signal from GRBs 
should be correlated in direction and time, greatly 
simplifying background rejection. All results re- 



ported here are based on data collected between 
April and October of 1997. Due to limitations in 
the data acquisition and archiving system at that 
early phase of operation, the livetime ranges be- 
tween 130 and 140 days, depending on the details 
of the analysis. 
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Figure 1. Survey of j/ + 17 flux predictions 

from cosmic accelerators of high energy neutrinos 
(adopted from the recent review by Learned and 
Mannheim poll. The atmospheric neutrino fluxes 
are from ref. |pQ|, for both vertical (1) and hori- 
zontal (2) fluxes within a circle deflned by a half 
angle of 4 degrees ( labeled "4 deg bin"). The 
curves do not include the normalization uncer- 
tainty, possibly 20% in magnitude. Numbered 
lines: (3) AMANDA-BIO hmit reported here; (4) 
Ref. |jll| for the core emission from 3C273 due to 
pp interactions. It also represents neutrino emis- 
sion from the AGN Blazar Mk501 during 1997 if 
it emits half of its TeV gamma ray flux in neu- 
trinos. (5) Crab Nebula. Model I from ref. |1^ , 
(6) Coma cluster according to ref. |1J], (7) core 
emission from 3C273 due to p-7 interactions [|l2[, 
(8)model||l^ for the relativistic jet of 3C273 in- 
cluding p-p and p-7 interactions, and Supernova 
remnant gamma-Cygni (9) and IC 444 (10) ac- 
cording to ref. |16| . The energy bounds on the 
AMANDA limit are restricted to the approximate 
region of sensitivity of the detector. 



2. ANALYSIS PROCEDURE - SEARCH 
FOR POINT SOURCES 

The various physics objectives are best ad- 
dressed by speciahzed analyses, although most 
are closely related. The search for sources of point 
emission can be used to illustrate the general 
methods. Trigger events are dominated by down- 
going atmospheric muons, so analysis techniques 
were developed to reject this background while 
retaining good efficiency for upgoing neutrino- 
induced muons. Unlike many neutrino detec- 
tors, the effective sensitivity depends strongly as 
a function of the background rejection require- 
ments, which are considerably weaker for searches 
for point sources than diffuse emission. 

The analysis procedure utilizes two essential 
characteristics of the signal to simplify the anal- 
ysis relative to atmospheric neutrino measure- 
ments. First, the sources are assumed to be point 
sources in the sky, so only events within a selected 
angular region are considered. Secondly, we use 
the topological characteristics of spectrally hard 
neutrino signal to reject poorly reconstructed 
atmospheric muons and atmospheric neutrinos, 
both of which have softer spectra. Topological 
variables include an estimate of muon energy and 
an assessment of the spatial fluctuation of the de- 
tected signals in a given event. The complete 
suite of variables was able to differentiate signal 
events from several classes of background topolo- 
gies. 

Monte Carlo based simulation programs deter- 
mined the effective area for background and neu- 
trino -induced muons. Several important results 
from these programs were tested by comparing 
the background simulation to the experimental 
data at various steps along the analysis chain. 

An iterative analysis procedure was developed 
to maximize the S/VBG, where the signal, S, was 
computed with an energy spectrum proportional 
to E~^ for the source. BG is background from 
atmospheric muons. After optimizing the analy- 
sis parameters, the sensitivity was evaluated for 
power law spectra with indices between 2.0 and 
3.0. 

The space angle resolution is determined from 
simulation. The upper panel of Figure shows 



that the median resolution is 3 degrees, and the 
lower panel indicates that this value only weakly 
depends on energy. Two studies were used to 
check the angular resolution and absolute off- 
set. First, events that simultaneously trigger 
the GASP ACT Hi and AMANDA provide a 
"test beam" containing single muons with direc- 
tional information provided by GASP. To im- 
prove the statistical accuracy of the investigation, 
a second study involved events which simultane- 
ously trigger the SPASE air shower array p6[ and 
AMANDA. Although the interpretation of these 
special events is complicated by the presence of 
multiple muons, which tend to reconstruct with 
worse angular precision than single muon events, 
the response of the detector to these events ap- 
pears to be correctly modeled. 



xlO 



^ 


0.7 


c 




^ 


06 










f. 


O.s 










u 


0.4 


■^ 






0.3 




0.2 




0.1 







'? 


6 




5 


n 






4 










^ 


3 


^ 




< 


2 



-• 

: -•- 


*** 

-•- 
-•- 
-•- 






' Mean ' 
RMS 




'4.024^ 
3.305^ 


Median 




3.078 J 


•- 

r(A) 


1 , , 


-• 
, 1 , 


-• 


-•- 


1,1,1 


^^^T^T" 


^^T»t«^r«^^» 



12 14 16 

A^ (Degrees) 



(B) 



2.5 



3.5 



4.5 



5 5.5 6 

log(EyGeV) 



Figure 2. Error in the space angle for simulted 
signal events with energy spectra proportional to 
i?i/^ . Top: distribution of error averaged over 
declination; Bottom: Space angle error as a func- 
tion of neutrino energy. 



The point source analysis yields an event sam- 
ple of 1097 events which are distributed on the 
sky as shown (Fig||): 
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Figure 3. Sky distribution of 1097 events in point 
source analysis. Coordinates are Right Ascension 
(RA) and declination (dec). 



Figure 4. The effective area for muon detection 
as a function of zenith angle for E^ between 0.1 
TeV and 100 TcV (180°is vertically up in local 
detector coordinates). 



Guided by the estimate of angular resolution, 
the sky was divided into 319 non-overlapping an- 
gular bins. The distribution of counts per sky 
bin is consistent with random fluctuations, which 
were determined by selecting all events within 
a declination band and randomly redistributing 
them in Right Ascension. 

The neutrino limits were computed according 
to 
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e// 



where A 



eff^ 



is the neutrino effective area 



weighted by the assumed neutrino energy spec- 
trum. This quantity is related to the muon effec- 
tive area shown in Figure |j. The factor T^ye is the 
livetime, and e is the efficiency due to finite an- 
gular resolution and also accounts for non-central 
source placement within an angular bin. The 
term ^{Nh, Nq) generates the 90% CL accord- 
ing to Feldman and Cousins p3] for signal events 



given the measured number of events in the bin, 
A^o and the expected background Nb determined 
from the events in the declination band contain- 
ing the source bin. The results of this calculation 
are shown in Figure ^ for various assumed spec- 
tral indices. 

The inferred limits on neutrino flux apply 
to point sources with continuous emission (or 
episodic emission averaged over a time interval 
of approximately 0.6 years) and power law en- 
ergy spectra with a fixed spectral index above the 
energy threshold of the detector. The limits for 
sources at large positive declination are compara- 
ble to the best published limits in the Southern 
skyi. 

3. ATMOSPHERIC NEUTRINOS 

The final event sample in the search for point 
sources contains both atmospheric neutrinos and 
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in the atmospheric neutrino sample from known 
physics effects is small (< 15%) for values of the 
event quality parameter greater than five. Back- 
ground simulations with much greater statistical 
precision are currently underway. 
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Figure 5. Preliminary neutrino flux limit (90% 
CL) on point sources of high energy neutrinos as 
a function of declination, averaged over RA. The 
limit is computed for a lower energy threshold is 
10 GeV. Note that the power law exponent refers 
to the neutrino energy spectrum. Also, neutrino 
absorption by the earth is taken into account. 



poorly reconstructed downgoing muons. The 
fraction of atmospheric v in the sample can be 
enhanced at the expense of sensitivity. Experi- 
mental data is dominated initially by background 
events - typically downward going atmospheric 
muons with poorly known directions. This can 
be seen in figure ^. as indicated by the flat be- 
havior for less restrictive selection criteria (qual- 
ity < 4.2). As selection criteria become progres- 
sively more restrictive (increasing values along 
the X-axis), the asymptotic flattening of the ra- 
tio (experimental data)/(Signal MC Atm.jy) indi- 
cates that the evolution of the experimental data 
becomes consistent with signal expectation in the 
vicinity of the plot where the (BG MC)/Exp ratio 
approaches zero. From this evidence (and visual 
inspection), we conclude that the contamination 
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Figure 6. Ratios of passing rates for sim- 

ulated background (BG MC), simulated atmo- 
spheric neutrino signal (Signal MC atm. v)^ and 
reconstructed experimental data (Exp) as a func- 
tion of " event quality" , a variable which measures 
the severity of the selection criteria. 



Figure shows that the angular distribution 
of 188 events is also consistent with the simu- 
lated distribution of atmospheric v events. Due to 
the elongated cylindrical geometry of AMANDA- 
BIO, the acceptance shows strong dependence on 
zenith angle. 

Thus, the angular dependence of the atmo- 
spheric neutrino sample is consistent with ex- 
pectation and contamination from background is 
small. Finally, the distribution of the number of 



OMs in an event is also consistent with expecta- 
tion (see figure H). 





- ' ' ' ' 


1 , , , 1 , , , 1 , , , 1 , , 


1 1 1 1 1 1 1 1 1 1 1 1 1 1 _ 


50 


^ 






• Data 


40 


- 






m^ AtmosMC 




H:,,^ 


V -i 


- 


30 


- 


m~ 


1 


-i- ; 


20 




.^ 


10 


- 


T T 


'^^ 1 - 








K>J»<1 


n 


- 





1 1 1 1 1 1 1 1 1 1 l^r l-r- 



-1 -0.9 -0.8 -0.7 -0.6 -0.5 -0.4 -0.3 -0.2 -0.1 

Cos(Zenith) 



Figure 7. Reconstructed zenith angle distribu- 
tion. The simulated atmospheric neutrino events 
(shaded boxes) are normalized to data (filled cir- 
cles). The vertical widths of the boxes indicate 
the errors computed using binomial statistics. 



4. WIMP SEARCH 

There is strong evidence for a non-baryonic 
component in the dark matter in the Uni- 
verse. Supersymmetric extensions to the Stan- 
dard Model provide promising dark matter can- 
didates, such as the neutralino. Dark matter in 
the galactic halo can interact with nuclei in the 
earth, lose energy and become trapped. The an- 
nihilation of neutralinos within the core of the 
earth will generate high energy neutrinos in the 
nearly vertical direction (cos(9) < —0.97). For 
this highly restricted search cone, the number of 
events observed agrees with prediction (17 pre- 
dicted, 15 observed). Therefore, the lack of a sta- 



tistically significant excess of neutrino events in 
the nearly vertical direction can be used to con- 
strain the flux of neutrinos from neutralino anni- 
hilation. Figure ^compares the AMANDA limits 
with predictions! 22 1 for a broad class of supersym- 
metric models, and with curves derived |23|] from 
published limits by MACRO § and BaksanQ, 
illustrating the potential of the technique. Sys- 
tematic errors associated with detector sensitiv- 
ity and neutrino oscillation physics are currently 
under investigation, but are expected to weaken 
the limits. 
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Figure 8. High energy muon flux predictions and 
experimental limits due to the annihilation of su- 
persymmetric particles in the center of the earth. 
See ref. |p3| for explanation of symbols. The 
AMANDA limit (90% C.L.) was corrected to cor- 
respond to a threshold of 1 GeV. 



5. HIGH ENERGY DIFFUSE FLUX 

Many models of diffuse emission from unre- 
solved sources predict an energy spectrum that is 
much harder than the atmospheric neutrino spec- 
trum. Therefore, a limit on the diffuse flux can be 
extracted by searching for an excess of events at 
large energies. A simple, although not very pre- 
cise, measure of energy of the neutrino-induced 
muons is counting the number of hit channels 
- optical modules that detected cherenkov light. 
Figure M shows the distribution of the hit mul- 
tiplicity for experimental data, for simulated at- 
mospheric neutrinos, and for an arbitrarily nor- 
malized, relatively hard energy spectrum. The 
experimental data agree with the atmospheric 
neutrino spectrum. From the non-observation of 
an excess of high multiplicity events, we derive 
an upper limit on an assumed diffuse E~^ spec- 
trum. The preliminary limit (90% CL) is of order 
dN/dE^ < IQ-^E^-^cm-'^s-^sr-^GeV-^. This 
experimental limit is below several models (e.g., 
ref. [0), although still above the recent bound 
established by Waxman and Bahcall[Q. The ex- 
perimental limit can be improved by developing 
better energy estimators using charge and topo- 
logical distributions. The response of the detec- 
tor to muons with energies in excess of 100 TeV 
is still under investigation. 
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Figure 9. Distribution of the number of optical 
modules for the data presented in Fig. 0. The 
MC simulation for atmospheric and a generic 
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6. GRB SEARCH 

Gamma Ray Bursts (GRBs) are thought to be 
produced by highly relativistic outflows originat- 
ing from a compact, explosive event. A search 
for emission of high energy neutrinos correlated 
in time and direction with GRBs was performed. 
We examined a total of 78 GRBs detected by 
BATSE. The Tgo duration of the bursts, defined 
as the time required to accumulate between 0.05- 
0.95 of the total counts, varied between 0.2s and 
200s. Guided by several models for neutrino emis- 
sion QQ, the analysis procedure optimized the 
effective area of the AMANDA detector for neu- 
trino energies between 100 TeV and 1 PeV. The 
lack of a statistically significant excess from any 
GRB examined in 1997 leads to a range of lim- 
its (the variation is primarily due to attenuation 



by the earth which depends on arrival direction) 
on the flucncc that span from Ej^dN/dE^, < 
0.05 - 1.0 X 10~i X mm(l, E^/Ebreak)TeVcm-^. 
The parameter Etreak, which is of order 500 TeV, 
is the energy that characterizes the change in 
the power law of neutrino spectrum (approxi- 
mately E~^ for El, < Ebreak and E~^ for neu- 
trinos with energies larger than Ebreak- Another 
change in the spectral index is expected at even 
higher energies, but this feature does not impact 
our analysis. Using all of the GRBs, a prelim- 
inary cumulative upper limit of E^, dN/dE^ < 
4x 10^''mm(l, E^ / Ei,reak)TeVcm'~'^ was also de- 
rived |l|l. For 
comparison, Waxman and Bahcali p7|| predict a 
fluence of 4.8 x 10"'^mm(l, Ei,/Ebreak)TeVcm~'^. 
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